Levy Distribution of Single Molecule Line Shape Cumulants in Low Temperature 

Glass 



O 

o 
o 

(N 

X) 
D 

m 



o 



X3 
O 

o 



> 
in 

o 

(N 
O 

o 
o 



I 

o 
o 



X 
S3 



E. Barkai", R. Silbey" and G. Zumofen'' 
" Department of Chemistry and Center for Material Science and Engineering, Massachusetts Institute of 

Technology, Cambridge, MA 02139. 
^ Physical Chemistry Laboratory, Swiss Federal Institute of Technology, ETH-Z, CII-8092 Zurich. 

(February 1, 2008) 



Abstract 

We investigate the distribution of single molecule line shape 
cumulants, ki, K2, ■ ■ •, in low temperature glasses based on the 
sudden jump, standard tunneling model. We find that the 
cumulants are described by Levy stable laws, thus generalized 
central limit theorem is applicable for this problem. 

Pacs: 05.40.-a, 05.40.Fb, 61.43.FS, 78.66.Jg 

Recent experimental advances Q have made it possi- 
ble to measure the spectral line shape of a single molecule 
(SM) embedded in a condensed phase. Because each 
molecule is in a unique static and dynamic environment, 
the line shapes of chemically identical SMs vary from 
molecule to molecule |^ . In this way, the dynamic prop- 
erties of the host are encoded in the distribution of sin- 
gle molecule spectral line shapes [0-^. We examine the 
statistical properties of the line shapes and show how 
these are related to the underlying microscopic dynami- 
cal events occurring in the condensed phase. 

We use the Geva-Skinner j|] model for the SM line 
shape in a low temperature glass based on the sudden 
jump picture of Kubo and Anderson |l^Jll|]. In this 
model, a random distribution of low-density (and non- 
interacting) dynamical defects [e.g., spins or two level 
systems (TLS)] interacts with the molecule via long range 
interaction (e.g., dipolar). We show that Levy statistics 
fully characterizes the properties of the SM spectral line 
both in the fast and slow modulation limits, while far 
from these limits Levy statistics describes the mean and 
variance of the line shape. We then compare our analyt- 
ical results, derived in the slow modulation limit, with 
results obtained from numerical simulation. The good 
agreement indicates that the slow modulation limit is 
correct for the parameter set relevant to experiment. 

Levy stable distributions serve as a natural general- 
ization of the normal Gaussian distribution. The impor- 
tance of the Gaussian in statistical physics stems from 
the central limit theorem. Levy stable laws are used 
when analyzing sums of the type '^Xi, with {xi} be- 
ing independent identically distributed random variables 
characterized by a diverging variance. In this case the or- 
dinary Gaussian central limit theorem must be replaced 
with the generalized central limit theorem. With this 
generalization. Levy stable probability densities, L-y^^(x), 



replace the Gaussian of the standard central limit theo- 
rem. Khintchine and Levy found that stable characteris- 



tic functions, Z/^^^(fc), are of the form |12 
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for < 7 < 2 (for the case r/ 7^ 0, 7 = 1 see |T^). Four 
parameters are needed for a full description of a stable 
law. The constant 7 is called the characteristic exponent, 
the parameter is a location parameter which is unim- 
portant in the present case, 2^ > is a scale parameter 
and — 1 < ?7 < 1 is the index of symmetry. When rj = Q 
[r] = ±1) the stable density L^^jiix) is symmetrical (one 
sided). Levy statistics is known to describe several long 
range interaction systems in diverse fields such as astron- 
omy [|l2 [^urbulence |l^ and spin glass 1 14 1 . Stoneham's 



theory 1 151 of inhomogeneous line broadening in defected 
crystal, is based on long-range forces and parts of it can 
be interpreted in terms of Levy stable laws Stone- 
ham's approach [|l^ is inherently static, while the SM 
line shape model considers both dynamical and static 
contributions from the defects. 

An important issue is the slow and fast modulation 
limits [|6|,pl|. Briefly, the fast (slow) modulation limit is 
valid if important contributions to the line shape are from 
TLSs which satisfy v ^ K {v K), where v is the fre- 
quency shift of the SM due to SM-TLS interaction and K 
is the transition rate of the flipping TLS (see details be- 
low). In the fast modulation limit, all (or most) lines are 
Lorentzian with a width that varies from one molecule to 
the other. For this case, the (Levy) distribution of line 
widths fully characterizes the statistical properties of the 
lines. The second, more complicated, case corresponds 
to the slow modulation limit. Then the SM line is typ- 
ically composed of several peaks (splitting) and is not 
described well by a Lorentzian. If a SM shows splitting, 
one can investigate the validity of the standard tunneling 
model of glass |l^ in a direct way, since the splitting of 
a line is directly associated with SM-TLS interaction . 
As mentioned, we demonstrate the existence of a slow 
modulation limit in SM-glass system. 

Following [|| we assume a SM coupled to non iden- 
tical independent TLSs at distances r in dimension d. 
Each TLS is characterized by its asymmetry variable 
A and tunneling element J. The energy of the TLS is 
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E = \J ^ J"^. The TLSs are coupled to phonons or 
other thermal excitations such that the state of the TLS 
changes with time. The state of the n th TLS is de- 
scribed by an occupation parameter, ^^(i), that is equal 
or 1 if the TLS is in its ground or excited state respec- 
tively. The probability for finding the TLS in its upper 
f = 1 state, p, is given by the standard Boltzmann form 
p = 1/{1 -|- exp[£^/(fcf,T)]}. The transitions between the 
ground and excited state are described by the up and 
down transition rates Ku^ Kd, which are related to each 
other by the standard detailed balance condition. 

The excitation of the n th TLS shifts the SM's tran- 
sition frequency by Vn- Thus, the SM's transition fre- 
quency is 



(2) 



where A^act is the number of active TLSs in the system 
(see details below) and ojq is the bare transition frequency 
that differs from one molecule to the other depending on 
the local static disorder. We consider a wide class of 
frequency perturbations 
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a is a coupling constant with units [Hz nm*] , 5* (fi) is 
a dimensionless function of order unity, is a vector of 
angles determined by the orientations of the TLS and 
molecule (in some simple cases depends on polar an- 
gles only), f{A, J) > is a dimensionless function of 
the internal degrees of freedom of the fluctuating TLS, S 
is the interaction exponent. The line shape of the SM is 
given by the complex Laplace transform of the relaxation 
function 
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provided that the natural life time of the SM excited 
state is long. The relaxation function of a single TLS 
was evaluated 1 11 1 based on methods developed in ||lffl 
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cosh (nt) + ^ sinh (fit) 



(5) 



with n ^ [X74 - Z/V4 -i(p- 1/2) i^i^T] 1/2^ S = f - 
i (p — ^) V and K = K^+Kd. For a bath of TLSs the hne 
shape, Eq. (^, is a formidable function of the random 
TLS parameters (r, 51, A, J) as well as the system param- 
eters (a,T, etc). In the fast modulation Hmit K :s> 
one finds a simpler behavior: all lines are Lorentzian with 
half width 



n 
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which varies from one molecule to the other. Eq. (g) 
shows the well known phenomena of motional narrow- 
ing. In the slow modulation limit K <^ \v\ one finds 
$(t) = 1 — p + pe~'^'^* implying that the line shape of 
a molecule coupled to a single TLS is composed of two 
delta peaks, the line shape of a molecule coupled to two 
TLSs is composed of four delta peaks, etc (splitting). 

The spectral line is characterized by its cumulants 
Kj (j — 1,2,---) that vary from one molecule to the 
other, and we investigate the cumulant probability den- 
sity P{Kj). We have derived the cumulants of the SM 
line shape, and the first four cumulants are presented in 
Table 1 [|l^ . We observe that cumulants of order j < 2 
are real while generally cumulants of order j > 2 are com- 
plex, implying that the moments of the line shape diverge 
when j > 2. The summation, X^n' Table 1 is over the 
active TLSs, namely those TLSs which flip on the time 
scale of observation r (i.e., Kn > 1/t). We consider the 
slow modulation limit, soon to be justified, which means 
that we consider the case Kn ^ I'n- To investigate this 
limit wc set Kn = in Table 1, then all the cumulants 
are real and are rewritten as kj — Hjni'-li, where Hj 
are functions of p only and Hi — p. Hi — p(\ — p), 
H3 = p(l — p){2p — 1) etc. Note that for k\ and K2 no 
approximation is made. 
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Table 1: Cumulants Kj of the SM line shape 



Let {■)rQAj denote an averaging over the random TLS 
parameters. The characteristic function of the j cumu- 
lant can be written in a form 



exp 
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PeS density of the active TLS and Bj ~ 

(27ra)^*J (fi) P {A, J) Hj. To derive Eq. (0) we have 
used the assumption of independent TLSs uniformly dis- 
tributed in the system. For odd j cumulants we find 



(exp [ikKj))^. 
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with characteristic exponent 7 — d/{5j) and the scale 
parameter 
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with Cj = cos (777/2) r(l — 7), ci = tt/2. Eq. (^) shows 
that odd cumulants are described by symmetrical Levy 
stable density, i.e., P{nj) = L~ffi[K,j). Two conditions 
must be satisfied for such a behavior, < 7 < 2 and 
/ sin [^-J (f^)] 0. The latter condition gives the sym- 
metry condition, 77 = 0, which means that negative and 
positive contributions to Kj are equally probable. 
For even cumulants and < 7 < 1 we find 

(exp {ikKj))^^j^j = L^,vi^) (10) 

with a scale parameter Eq. (^) and with Levy index of 
symmetry 

^ {PHA,J)\H,\-<)^^j ■ ^ > 

Eq. ( |To| ) implies that even cumulants are distributed 
according to P{Kj) = L^^,,(kj). We note that the asym- 
metrical Levy functions, with j] 7^ ±1,0, only rarely find 
their applications in the literature. The characteristic 
exponent 7 depends only on the general features of the 
model (namely on d and d). In contrast the Levy in- 
dex of symmetry 77 depends on the details of the model 
and on system parameters (T, etc). For j = 2 we have 
Hj — \Hj\ and then r/ = 1 so the Levy density is one 
sided, as is expected since K2 > 0. 

As mentioned, in the fast modulation limit, the ran- 
dom line width in Eq. (||) characterizes the statistical 
properties of the spectral lines. Using the approach in 
Eqs. (^[^) one can show that P{T) = ici/l25).i(f ) with 
the scale parameter 2^/(25) given by Eq. (|sj) with j = 2 
and H2 =p{l-p)/K. 

In what follows we exhibit our results and compare 
to simulations based on the standard tunneling model 
of low temperature glass [p^ . We use system parame- 
ters given by Geva and Skinner |^ to model terrylene in 
polystryrene. The SM-TLS interaction is dipolar, hence 
(5 = 3, and we consider spatial dimension c? = 3. The dis- 
tribution of the asymmetry parameter and tunneling ele- 
ment is P(A)P(J) = N-^J-^ for2.8xlO-^K < J < 18K 
and < ^ < 17K, N denoting a normalization con- 
stant. We use f{A, J) = A/E and define a TLS to be 
active ii K > l/r, r = 120 sec is the time of experi- 
mental observation. In this way the averaging (• • 
becomes r independent. The rate of the TLS is given 
by K ^ cJ^E coth{PE„/2) and c = 3.9 x lO^K-^Hz is 
the TLS phonon coupling constant. Additional system 
parameters are the coupling constant a = 3.75 x 10^^ 
nm^Hz and the TLS density 1.15 x 10^^ nm^^. Accord- 
ing to Eqs. (|^)-(|l^), only the scale parameter Zj depends 
on the orientation of the TLS and SM, through *(fi). It 



is therefore reasonable to assume simple forms for ^'(fi). 
We consider two examples, model 1 (Ml) for which '^{^l) 
is replaced with a two state variable (i.e., a spin model) 
5' = lor\E' = — 1 with equal probabilities of occurrence 
and model 2 (M2) ^'(f7) = cos(6l), with 9, the standard 
polar coordinate, distributed uniformly. With these def- 
initions we calculate the symmetry index rj and the scal- 
ing parameter z-y and compare between the theory and 
numerical simulation. 



0.35 - 




-10 -5 5 10 

K,/Z, 



FIG. 1. Scaled probability density of first cumulant 
P{fti)zi versus k\I Z\. Symbols are the simulation results ob- 
tained from 4000 molecules for different cases indicated in 
the figure. The theory, plotted as a solid curve, predicts a 
Lorentzian density P{ni) ~ Li,o(«:i) with a scaling parame- 
ter zi which varies from one set of data to the other. 

We consider the first two cumulants ki and K2 (i.e, 
the line shape mean and variance). Since d — S we 
find P{ki) = Lifi(Ki) which is the Lorentzian density, 
and P{n2) = ^1/2, 1(^*2) which is Smirnov's density. We 
have considered two temperatures for the two models Ml 
and M2. As shown in Fig. and ||, a scaling behavior 
is observed and all data collapse on the Levy densities 
Li,q{ki) and ii/2,i(K2) respectively. In Fig. |^ and || we 
have rejected TLSs within a sphere of radius rmin — 1 nm, 
demonstrating that our results are not sensitive to a short 
cutoff. Also shown in the inset of Fig. || is P(Re[K3]) 
which is distributed according to L1/3 o(Re[K3]) and a 
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scale parameter zi/3 given in Eq. (^. The Levy behav- 
ior of Ki, K2 and Re[K3] holds generally and is not limited 
to the slow modulation limit since these random variables 
do not depend explicitly on the rates K. 



so the assumption of slow modulation limit is justified. 
For model M2, we see slightly larger deviations between 
the theory and numerical results, due to the angular de- 
pendence of model M2, "^{Vl) = cos{6), which reduces 
the typical frequency shift \v\ compared to model Ml. 
We conclude that the present theory can be used as a 
criterion for the validity of the slow modulation limit. 
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FIG. 2. Same as Fig. (nl) for the second cumulant. We 
show P(k2)2i/2 versus K2/z1^2- The solid curve is Smirnov's 
density ii/2,i('^2). In the inset we show the same as Fig. (^) 
for Re[K3]], the solid curve is Levy density Li/3_o(Re[Ki/3]). 

Consider the distribution of Re[K4], which in the 
slow modulation limit is distributed according to 
Li/4_^(Re[K4]), Eq. (prj). The question remains if such a 
slow modulation limit is valid for the standard tunneling 
model parameters we are considering. The slow modula- 
tion limit is expected to work when K <C |!^|. For large 
enough r this inequality will fail; however, depending on 
system parameters, we expect that contributions from 
TLS situated far from the SM are negligible. We also 
note that according to the standard tunneling model the 
rates K are distributed over a broad range, albeit with 
finite cutoffs that insure that the averaged rate is finite. 
To check if the slow modulation limit is compatible with 
the standard tunneling model approach, we compare our 
slow modulation results with those obtained by simula- 
tion in Fig. ||. We also show simulation results in which 
all rates are set to zero {K = 0). For model Ml, we find 
that the deviation between simulation and theory is small 



FIG. 3. Same as Fig. (|l|) for Re[K4\ and for temperature 
T — 1.7K. The symbols are the simulation results obtained 
for four different cases as indicated in the figure and for 12000 
molecules. The solid curve is Li/4_^[Re(K4)] with an index of 
asymmetry rj = 0.6104. In the inset we show the power law 
tail of the scaled probability density on a log-log plot. 

Depending on system parameters, Levy statistics may 
become sensitive to the finite cutoff r^^^, Physically, the 
cutoff can be important since the power low interaction 
is not supposed to work well for short distances |^ . Our 
results were derived for rmin = 0, while for small though 
finite rmin one can find intermittency behavior, i.e., the 
ratio (k2)/(k2)^ (as well as similar dimensionless ratios) 
is very large. When rmin is large one finds a Gaussian 
behavior. The phenomena of intermittency in the con- 
text of a reaction of a SM in a random environment was 
investigated in |Q. Generally high order cumulants are 
more sensitive to finite cutoff and for results in Fig (||) 
Tmin = was chosen to see the proper decay laws in the 
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wing. 

To conclude, we showed that the generaUzed central 
limit theorem can be used to analyze distribution of cu- 
mulants of SM line shapes in glass. We note that besides 
cumulants, Levy statistics can be used to analyze other 
statistical properties of SMs in disordered media Q . 
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